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Lyman "bump" galaxies - I. Spectral energy distribution of 
galaxies with an escape of nebular Lyman continuum 
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ABSTRACT 

It is essential to know galactic emissivity and spectrum of Lyman continuum (LyC) 
in order to understand the cosmic reionization. Here we consider an escape of nebular 
LyC from galaxies and examine the consequent spectral energy distribution. It is usu- 
ally assumed that hydrogen nebular LyC mostly produced by bound-free transitions 
is consumed within photo- ionized nebulae (so-called "on-the-spot" approximation). 
However, an escape of the continuum should be taken into account if stellar LyC es- 
capes from galaxies through "matter-bounded" nebulae. We show that the escaping 
hydrogen bound-free LyC makes a strong bump just below the Lyman limit. Such a 
galaxy would be observed as a Lyman "bump" galaxy. This bump results from the 
radiation energy re-distribution of stellar LyC by nebulae. The strength of the bump 
depends on electron temperature in nebulae, escape fraction of stellar and nebular 
LyC, hardness of stellar LyC (i.e. metallicity, initial mass function, age and star for- 
mation history), and IGM attenuation. We can use the bump to find very young (~ 1 
Myr), massive (~ 100 M@), and extremely metal-poor (or metal- free) stellar popula- 
tions at z < 4. Because of the bump, 900 A-to-1500 A luminosity density ratio (per 
Hz) becomes maximum (2-3 times larger than the stellar intrinsic ratio) when about 
40% of the stellar LyC is absorbed by nebulae. The total number of escaping LyC 
photons increases due to the escape of nebular LyC but does not exceed the stellar 
intrinsic one. The radiation energy re-distribution by nebulae reduces the mean energy 
of escaping LyC only by » 10% relative to that of stellar LyC. Therefore, the effect 
of the escape of nebular LyC on the reionization process may be small. 

Key words: cosmology: theory — cosmology: observations — galaxies: evolution — 
galaxies: high-rcdshift — H II regions — intergalactic medium 



1 INTRODUCTION 

Understanding the reionization of the intergalactic medium 
(IGM) is one of the most important issues in cosmol- 
ogy. The end of the hydrogen reionization epoch is prob- 

ably at z « 6 which is inferred from QSOs (e.g., 

iBecker et all l200ll: iFan. Carilli. Keating! |2006ft . Lyman a 



emitter s (LAEs) (Kashikawa et al. 120061) . and Gamma-ray 



bursts l|Totani et al.1 120061 ) at the redshift. Future 21 cm 
tomography wi ll reveal the ionization history of hydrogen 
in detail (e.g.. iMadau. Meiksin. fc Reesl Il997l ). Even after 
we know the history, however, the problem of the ionising 
source may remain st ill open. We have not found many ob- 
jects at z > 6 so f ar |lve et alj|2006l ; ISalvaterra et al.ll2009l: 
iTanvir et afll2009h . Moreover, we do not know the Lyman 
continuum (LyC) emissivity of high-z objects. 

Galaxies are thought to be more probable source 
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for the hydrogen reioniz a tion than QSOs (e.g., 
IMadau. Haardt. fc Reed Il999l ). Ilnoue. Iwata. fc Deharvengl 
found that the LyC emissivity relative to non-ionising 
ultraviolet (UV) (or escape fraction of LyC) should increase 
by an order of magnitude from z < 1 to z > 4 so as to match 
the observed UV luminosity density of galaxies with the 
hydrogen ionization rate estimated from Lyman a forest. 
This trend is consistent with a remarkable contrast between 
z < 1 and z ~ 3 in direct observations of LyC from galaxies 



there is no detection of LyC at z < 1 dLeitherer et all 
Hurwitz, Jelinskv, fc Dixon! 1 19971: iDeharveng et al.l 



Siana et al.l 120071 : ICowie. Barger. fc Trouille 

" ' 20061 : 



1995: 



2001: 



2009), but 



Grimes et al.l 



a ma rginal exception ( Be rgvall et al.l 
2007), whereas there are direct detections of LyC from 
Lyman break galaxies (LBGs) an d LAEs at z ~ 3 
dSteidel. Pettini.'fc Adelbergerl 1200 ll : IShaplev et all 120061 : 



llwata et al. 20091 ). This evolving escape fraction is also 



found in a simulation of galaxy formation with the radiative 
transfer l|Razoumov fc Sommer- Larsenll2006l . l2007l ). 
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The latest observations with the Subaru te lescope have 
found LAEs emitting extremely strong LyC (jlwata et al.l 
2009). Indeed, some of them are brighter in LyC (~ 900 
A) than in UV (~ 1500 A). Such extreme "blue" colours 
cannot be explained by stand ard stellar populat ions with a 
Salpeter initial mass function (|lwata et al.ll2009l ). If we esti- 
mated the escape fraction for them, we would obtain a value 
larger than unity. This shows our poor knowledge of LyC 
emissivity and spectrum of galaxies. We should consider the 
LyC spectrum of star-forming galaxies more carefully. 

This paper examines effects of nebular emis- 
sion on the spectral energy distribution (SED) 
of galaxies. Such studie s have been made s o far 
by several authors (|Stasmska fc Leitherer 19961; 



Table 1. Models of stellar population and SED. 



Mov 



Rocca-Volmerange. Fiod l200ll; [ Zaekrisson et all 
200ll; ISchaererl l2002l; I Zackrisson. Bergvall. fc Leited 120081 ; 
Schaerer fc de Barroa [2009J). However, we take into ac- 



count, for the first time, an escape of nebular LyC from 
photo-ionized nebulae and galaxies. The nebular LyC is 
usually assumed to be absorbed wit hin the nebulae, the so- 
called "on-the-spot" approximation (lOsterbrock fc Ferland I 
2006). This is a very good treatment for "photon-bounded" 
nebulae from which no LyC escapes. If stellar LyC escapes 
from "matter-bounded" nebulae and escapes from galaxies 
as observed at z ~ 3, however, nebular LyC may also 
escape. We discuss this latter case here. 

The rest of this paper consists of four sections. In Sec- 
tion 2, we describe modelling procedures: stellar population 
models, calculations of nebular emission, and IGM attenu- 
ation. The predicted spectral energy distributions, 900 A- 
to-1500 A luminosity density ratios, and other results are 
presented in Section 3. In Section 4, we discuss validity of 
assumptions made in modelling and a few implications from 
the results. The final section is devoted to the conclusion. 



2 MODEL 

2.1 Stellar emission 

We use population synthesis models taken from literature. 
Parameters characterising SEDs of galaxies are metallicity 
and initial mass function (IMF) as well as star formation 
history and age. Moti vated by the disco very of LAEs emit- 
ting very strong LyC l|lwata et al.| [2009) , we consider stellar 
populations emitting LyC as strongly as possible. Thus, we 
consider very young age of 1 Myr after an instantaneous 
burst, except for §3.2 where we also consider a continuous 
star formation with a longer age. The metallicity assumed 
here is Z — 0.0004 which is 1/50 of the Solar metallicity 
and the lowest one available in the code Starburst99 ver- 
sion 5.1 l|Leitherer et al.ll"l 999) 3 Then, we consider two IMF 
slopes of —2.35 IjSalpeterl Il955f ) and —0.1 (extremely top- 
heavy). We call this two models A and B, respectively. The 
mass range is assumed to be 1-100 M© for both models. 
In addition, we consider two mo re cases: Z = 10 (ex - 
tremely metal-poor [EMP] stars; iBeers fc Christliebl 120051 ) 



1 The stellar track assumed is a Padova track with AGB stars. 
The atmosphere model assumed is "PAULDRACH/HILLIER" as 
the recommendation of the code. This choic e of the atmosp here 
model may affect the stellar LyC emissivity dSchaeredl2003l) . 
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Figure 1. Schematic picture of clumpy ISM. Dense regions which 
remain neutral against the ionization by the stellar radiation are 
shown by the thick shades and ionized nebulae formed at the 
surface of the dense regions are shown by the thin shades. Other 
space is filled with diffuse ionized gas which has negligible opacity 
for Lyman continuum. 



and Z = (metal-free or Population-Ill stars) with the 
Salp eter IMF s lope a nd 50-500 Mq . These models are taken 
from lSchaereH (2003) and called C and D, respectively. Table 
1 is a summary of the stellar population models considered 
in this paper. 

We have to note that modelling of stellar LyC is still 
developing. For example, stellar rotation, which is not taken 
into account here, may increase the LyC emissivity of stars 
with an age betwee n 3-10 Myr by a factor of about 2 
l|Vazquez et al.l 2007). Thus, stellar LyC models have such 
an amount of uncertainty. 



2.2 Nebular emission 

2.2.1 Structure of the ISM 

The interstellar medium (ISM) in galaxies is clumpy. There 
should be dense clumps which remain neutral against the 
ionising radiation from stars. Considering a typical density 
of neutral hydrogen (~ 1 3 cm" 3 ) and a typical size (~ 5 
pc) found in dark clouds (|Mverslll978h . we find 10 5 as the 
optical depth at the Lyman limit and 10 even at 100 A. 
Thus, all the LyC photons along a ray from an ionising star 
are absorbed if the ray intersects a dense clump. If the inter- 
clump medium is diffuse and ionized highly enough to have a 
negligible opacity for LyC0, the escape fraction of the pho- 
tons is determined by the covering fraction of the clumps 

2 If not, we could not observe strong LyC from galaxies. In ad- 
dition, shock ionization by multiple supernovae would contribute 
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around an ionising star. Thus, it becomes independent of 
the wavelength. Ionized regions are formed at the surface 
of the neutral clumps. The nebular emission from the re- 
gions includes LyC. The escape fraction of the nebular LyC 
is assumed to be the same as that of the stellar LyC for 
simplicity. This may be good in average sense when stars 
and clouds are well mixed. Figure 1 is a schematic picture 
of such a clumpy medium. 



2.2.2 Luminosity density 

Let us denote the escape fraction of LyC as / csc . This is the 
number fraction of LyC photons which escape from a galaxy 
relative to the photons produced in the galaxy. As described 
above, we omit the frequency dependence of / eS c- However, 
we will see the effect in Section 4.1. 

The ionization equilibrium of hydrogen in photo-ionized 
nebulae with / csc is 



Q„(l - /csc) + Qneb(l - /csc) = / n p n c a A (T c )dV 



(1) 



where Q* and Q nc b are the production rates of LyC pho- 
tons by stars and by nebulae, respectively, n p and n c are 
proton and electron number densities, respectively, and ola 
is the recombination rate to all states. The integral is per- 
formed over the entire volume of the nebulae, V. The re- 
combination rate depend s on the local electron tempera- 
ture in the nebulae, T B (|Osterbrock fc Ferlandl I2006T ). If 
we assume T c to be uniform in the nebulae, we can write 
J n p n B ctAdV = ola J n p n e dV. The production rate of neb- 
ular LyC photons can be expressed as Q n eb = a% J n p n c dV 
(for a uniform T c ) because LyC photons are produced by the 
recombination to the ground state (its rate is ai), except for 
the free- free contribution which is negligible in photo-ionized 
nebulae. Therefore, equation (1) is reduced to 



Q*(l-/csc) 

n p n c dV = ; — -, — { — 

P Qfl(Tc) +Ql(Tc)/ C! 



(2) 



Note that the Case B recomb ination rate is qb = oa — on 
jOsterbrock fc Ferland Il2006h . 

The luminosity density of the nebular emission at fre- 
quency v is 



Lncb 
u — 



7° eb (T e )n p n e dV . 



(3) 



where j™ ch is the nebular emission coefficient which does 
not only depend on the temperature T e but also on the den- 
sity. However, we can omit the density dependence here (see 
§2.2.3). From equations (2) and (3), we obtain 



neb/rj-, ■ 



.(1-/. 



a s (Tc) +ai(T e )/ 0! 



(4) 



for a uniform T c . Note that equation (4) is equivalent to 
the standard t reatment of neb ular emission, for example, 
equation (2) in Schacrcr (2002), if / CS c = 0. The evaluation 
of equation (4) needs recombination rat es. We obtain the 
rates from an interpolation of Table 3 in iFerlandl l|l980t) . 



to keep low neut r al fra ction in the inter-clump medium as shown 
bv lYaiima et al.1 j2009t) . 
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Figure 2. Hydrogen nebular emission coefficients for (a) electron 
temperature T e = 10 4 K, (b) 10 5 K, and (c) 10 6 K. The bound- 
free (dashed), free-free (dot-dashed), and two-photon (dotted) co- 
efficients are shown. The solid line in each panel is the sum of the 
bound-free and the free-free and it should be compared with the 
calculations by Ferland (1980) which are shown by small filled 
circles. 



2.2.3 Emission coefficients 

The nebular emission coefficient is the sum of three emission 
processes of bound-free, free-free, and two-photon as 

neb / 



7" eb (Tc) = T bf (Te) + 7 5 (To) + 7 ' q (Te) . 



(5) 



We consider, for simplicity, nebulae composed of hydro- 
gen only since the effect of helium is small as shown later 
(§4.2). The free-free and bound-free emission coefficients 
are c alculated by equations (3) and (9) of iMewe et al.l 
19861). The two-photon emission coefficient is calculated as 



Osterbrock fc Ferland Il200el ) 



„ 2 q 



2hyP{y)X 2 c l aB , 



(6) 



where h is the Planck constant, y is the frequency nor- 
malized by that of Lya, P(y) is the normalized spectral 
shape whose approximated fu nctional form is given by 
iNussbaumer fc Schmutzl (1 19841 ). and X 2q is the probabil- 
ity of being in 2 2 S state after one Case B recombin ation. 
We adopt X 2q = 0.32 jSpit zer fc Greensteinl Il95ll ). The 
two-photon emissivity depends on density if nebulae have a 
high density, but only the sm all density limit (< 10 4 cm" 3 ; 
lOsterbrock fc Ferland I [2006) would be enough for this pa- 
per. We also note that the Case B approximation, which 
is suitable for nebulae optically thick for Lyman series lines 
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Rest wavelength [ftm] 

Figure 3. Spectral energy distributions (SEDs) with escaping nebular Lyman continuum: (a)— (d) for models A— D in Table 1. The 
assumed parameters are shown in panel (a): electron temperature T e = 10 4 K in nebulae, escape fraction / osc = 0.5, and 1 Myr after an 
instantaneous burst. The stellar (dashed), nebular (dot-dashed), escaping stellar+nebular (solid), and hot gas (dotted; T c = 10 6 K and 
5% of the bolometric luminosity) SEDs are shown in each panel. The vertical axis is normalized by the mass of the star cluster. 



l|Osterbrock fc Ferland f2006l ) , does not conflict with our as- 
sumption of nebulae optically thin for LyC because the ab- 
sorption cross sections for Lyman lines are a few orders of 
magnitude larger than that for LyC. 

Figure 2 shows the nebular emission coefficients for 
T c = 10 4 , 10 5 , and 10 6 K: bound-free (dashed lines), free- 
free (dot-dashed lines), and two- photon (dotted lines). We 
have confirmed that sums of the bound-free and free-free 
coefficient s (solid l ines) agree with the values presented in 
Table 1 of lFerlandl (|l980l ) for A > 912 A (small filled circles) 
within about 10% difference. 



2.3 IGM attenuation 

Radiation with a wavelength shorter than Lya in the rest- 
frame of the source is absorbed by neutral hydrogen re- 
mained in the IGM. We use a Monte Carlo simulation of the 
IGM attenuation bv llnoue fc Iwatal l|2008l ). This simulation 
is based on an empirical distribution function of the column 
density, the number density, and the Doppler parameter of 
absorbers in the IGM which is derived from the latest ob- 
servational statistics. The predicted Lya decrements agree 
with the observations in the full range of z = 0-6 excellently. 



2.2. 4 Hot gas heated by supernovae and stellar winds 

In addition to photo-ionized nebulae, shock-ionized nebu- 
lae produced by multiple supernovae and st ellar winds may 
contri bute to LyC emission. According to iLeitherer et al.l 
(1999), the mechanical luminosity produced by supernovae 
and stellar winds is ~ 1% of the bolometric luminosity of 
stars and well less than 10% of it. We assume here 5%, a 
somewhat large value, to show that this component is less 
important. We also assume the temperature to be 10 6 K for 
this component. 



3 RESULT 

3.1 Spectral energy distribution 

First, we show the resultant SEDs without IGM attenua- 
tion in Figure 3. Panels (a)-(d) correspond to models A- 
D in Table 1. The assumed quantities are T c — 10 4 K, 
/esc = 0.5, and age of 1 Myr after an instantaneous burst. 
The dashed lines are the stellar SEDs. We can see that 
LyC becomes harder from the model A to D (see also Ta- 
ble 3) . The dot-dashed lines are nebular emissions produced 
in photo-ionized nebulae. We can see bound-free Lyman, 
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Figure 4. Close-up spectra around the Lyman limit: (a)-(d) for 
models A-D in Table 1. The assumed parameters are shown in 
panel (c): electron temperature T c = 10 4 K in nebulae, escape 
fraction / csc = 0.5, and 1 Myr after an instantaneous burst. The 
stellar (dashed) and escaping stellar+nebular (solid) spectra are 
the same as Figure 3. The dotted lines show spectra affected by 
IGM mean attenuation. Redshift of the galaxy is z = 1 to 6 from 
top to bottom in each panel. The vertical axis is normalized by 
the flux density at Lya. 



Balmer, and Paschen continua and two-photon continuum. 
Free-free emission is also taken into account but negligible 
for this case. The solid lines are escaping stellar+nebular 
emissions: sum of the two emissions but reduced by a factor 
of /esc in LyC. We can see a significant contribution of the 
nebular emission. In particular, the bound-free LyC makes a 
bump just below the Lyman limit for models B-D. Even for 
the model A, the escaping LyC is a factor of 2 larger than 
the stellar one at ~ 900 A. 

The significant contribution of the nebular LyC is due 
to the energy re-distribution by nebulae. When a LyC pho- 
ton ionizes a hydrogen atom, the photon energy moves to 
a photo-electron. The energy is thermalized in nebulae col- 
lisionally. Then, an electron recombines with a proton. If 
the recombination is to the ground state, we get a photon 
with an energy of 13.6 eV + kinetic energy of the electron. If 
T c ~ 10 4 K, the kinetic energy is an order of 1 eV. Thus, the 
photon wavelength is not far from the Lyman limit. More- 
over, the recombination probability increases for an electron 
with a lower kinetic energy. Therefore, nebular LyC shows 
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Figure 5. Effect of electron temperature T c on the escaping 
SEDs: (a) for the model A and (b) for the model D. The assumed 
escape fraction is / osc = 0.5 and the age is 1 Myr after an instan- 
taneous burst. Four cases are shown in each panel: T c = 5 X 10 3 
K (dot-dashed), 1 X 10 4 K (solid), 2 X 10 4 K (dashed), and 5 X 10 4 
K (dotted). 



a peak at the Lyman limit. This is already found in the 
emission coefficient shown in Figure 2. 

The dotted lines in Figure 3 are emissions from hot gas 
produced by supernovae and stellar winds. We can conclude 
that this componen t is less important On the other hand, 
Yaj ima et all l|2009h have shown an importance of shock ion- 
ization for a large escape fraction; shock ionization in addi- 
tion to photo-ionization keeps the content of neutral hydro- 
gen lower and suppresses the LyC opacity through the ISM. 
This is not conflict with our argument that LyC emitted by 
shocked gas itself is negligible. 

Figure 4 shows model spectra affected by IGM attenu- 
ation (dotted lines). The assumed parameters are the same 
as Figure 3, and thus, the stellar (dashed) and escaping stel- 
lar+nebular spectra (solid) are the same as Figure 3. The 
spectra are normalized by the flux density at Lya and the 
unit of the vertical axis is arbitrary if it is for flux density 
(per Hz). We can see a prominent "bump" just below the 
Lyman limit due to the nebular bound- free LyC. Neutral hy- 
drogen remained in the IGM absorbs radiation below Lya . 
We apply IGM mean attenuation bv llnoue fc Iwatal (|2008l ) 
as described in §2.3.1. Each dotted line corresponds to dif- 
ferent redshift at which the object lies: z = 1, 2, 6 from 
top to bottom. We find that the Lyman limit bump is visi- 
ble up to z — 3 in this figure. Spectroscopic or narrowband 
observations tracing this wavelength range for z < 3 is quite 
interesting to confirm the reality of the model. If the object 
lies along a sight-line less opaque than the average, we may 
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Figure 6. Luminosity density ratio (per Hz) of Lyman contin- 
uum (900 A) to non-ionizing ultraviolet (1500 A) as a function of 
escape fraction. The assumed parameters are electron tempera- 
ture T e = 10 4 K and 1 Myr after an instantaneous burst. Models 
A— D in Table 1 are shown with each line style. 

find the bump from z 4 but not from z > 5 (see Figs. 8 
and 11 in llnoue fc Iwatall2008t ). 

Figure 5 shows the effect of the assumed electron tem- 
perature T c on the escaping stellar+nebular SEDs: T c = 
5 x 10 3 K (dot-dashed lines), 1 x 10 4 K (solid), 2 x 10 4 K 
(dashed), and 5 x 10 4 K (dotted). Other parameters (/esc 
and age) are the same as Figure 3. We show the models A 
and D only but the models B and C are qualitatively similar. 
We find that the shape of the Lyman limit bump becomes 
more peaked for lower T c . This is because the number of 
electrons with larger kinetic energy is smaller for lower T c in 
the Maxwell-Boltzmann distribution, thus, the photon en- 
ergy emitted by recombination becomes closer to the Lyman 
limit for lower T e . The same thing is true for the behavior 
around the Balmer limit. 

3.2 900 A-to-1500 A luminosity density ratio 

When estimating the escape fraction of LyC from di- 
rect observations, we often have to assume an in- 
trinsic stellar ratio of LyC to non-ionizing U V lu- 
minosity densities (see §4 of | lnoue et all 120051 ) For 
example, ISteidel. Pettini. fc Adelbergerl ( 200lh assumed 
Lu90o/ Luiboo = 0.33 for their z ~ 3 LBGs and ISiana et al.l 
(2007) argued that 

Lvqoo I -ki/isoo = 0.17 is better for their 
Z ~ 1 galaxies based on SED fitting. However, escaping stel- 
lar+nebular SEDs shown in Figure 3 are completely differ- 
ent from the stellar intrinsic ones because of the Lyman limit 
bump by the nebular bound-free LyC. Here, we present the 
900 A-to-1500 A luminosity density ratio of escaping stel- 
lar+nebular spectra. 

Figure 6 shows the 900 A-to-1500 A luminosity den- 
sity ratio as a function of escape fraction / osc for the case 
with T c = 10 4 K and 1 Myr after an instantaneous burst. 
Four models in Table 1 are shown by each line style: solid 

3 In literature, 1500 A-to-900 A ratio is more popular, but we 
here take the inverse to avoid the divergence when / osc — » 0. 



Figure 7. 900 A-to-1500 A luminosity density ratio (per Hz) as a 
function of the duration of constant star-formation for models A 
and B. For model A, stellar ratio is the dashed line and escaping 
stellar+nebular ratio is the solid line. For model B, stellar ratio 
is the dotted line and escaping stellar+nebular ratio is the dot- 
dashed line. The electron temperature T c = 10 4 K and escape 
fraction / osc = 0.5 are assumed. Two gray thick marks indicate 
the values assumed in Steidel et al. (2001, SPA01) and Siana et 
al. (2007, S+07). 

for model A, dashed for model B, dot-dashed for model C, 
and dotted for model D. At / csc = 1, there is no nebula, 
and we see the stellar intrinsic ratio for each model. As / eS c 
decreases towards 0, the ratio increases, presents a peak, de- 
creases, and reaches 0. The peak appears at / esc = 0.58-0.62 
and the peak ratio is a factor of 1.7-2.8 larger than the stel- 
lar ratio. Note that the 900 A-to-1500 A luminosity density 
ratio becomes maximum not when 100% of the stellar LyC 
escapes but when 40% of it is absorbed by nebulae. This 
is the Lyman bump effect due to the nebular bound-free 
emission. 

Figure 7 shows the effect of the age on 900 A-to-1500 A 
luminosity density ratio. We again assume T e = 10 4 K and 
/csc = 0.5 and also assume a constant star formation for this 
analysis. We only show the results of the models A and B. 
Overall feature found in Figure 7 is that the ratio decreases 
after the first a few Myr and reaches an asymptotic value. 
This is because 900 A radiation comes only from very mas- 
sive stars whose life-time is less than a few Myr but 1500 A 
radiation comes from not only such massive stars but also 
relatively long-lived intermediate mass stars. Thus, the 1500 
A luminosity density increases until the intermediate mass 
stars have turnover. The contribution of such intermediate 
mass stars to 1500 A radiation depends on the IMF. For the 
model B, the contribution is small, and then, the luminosity 
density ratio reaches the asymptotic value earlier than the 
model A. For models C and D, the ratio is almost indepen- 
dent of the age because there is no intermediate mass stars 
in these models. 

Let us compare the luminosity density ratio predicted 
here with those assumed in literature which are shown 
by gray thick marks in Figure 7. T he ratio proposed 
bv lSteidel. Pettini. fc Adelbergerl (|200ll ) corresponds to the 
stellar intrinsic ratio (dashed line) for 10-20 Myr age of the 
model A but to the escaping stellar+nebular ratio (solid) for 
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— T e =5 10 J K 
1 10 4 K — 

3 10* K 

5 10 4 K — 



J I I I I I I I I I I I I I I I I I i_ 



Table 2. Coefficients for /3 in equation (11) and mean energy of 
nebular Lyman continuum. 



0.0 0.2 0.4 0.6 0.8 1.0 
Escape fraction 

Figure 8. Ratio of the stellar+nebular escape rate to the stellar 
production rate of LyC photons, or "effective" escape fraction, as 
a function of escape fraction. Four lines show different electron 
temperature T : 5 X 10 s K (dot-dashed), 1 X 10 4 K (solid), 2 X 10 4 
K (dashed), and 5 X 10 4 K (dotted). 



100 M yr age of the model. The ratio proposed bv lSiana et alJ 
corresponds to the stellar intrinsic ratio for > 100 
Myr age of the model A. On the other hand, we expect 
much larger ratios for younger age: I/„9oo/Li;i500 = 1-1 
for the model A with the nebular contribution if age less 
than a few Myr. This ratio i s a f actor of 3 larger than 
ISteidel, Pettini, fc Adelbergerl <|200ll ) and a factor of 6 larger 
than lSiana et alJ (|2007l ). If we change the IMF (model B) or 
metallicity and mass range (models C and D), even larger 
ratios are expected (see Figure 6). 

3.3 Total number of escaping LyC photons 

We have seen that the 900 A luminosity density increases 
due to the nebular bound-free emission. This is the radia- 
tion energy re-distribution by nebulae. How about the total 
number of photons in escaping LyC? The LyC photon es- 
cape rate is given by 



Qcsc = (Q* + < 



(7) 



with the escape fraction of stellar and nebular LyC, /esc, and 
the production rates of LyC photons by stars, Q*, and by 
nebulae, Q nc b. From the ionization equilibrium in equation 
(1) combined with equation (2), equation (7) is reduced to 

Qcsc CVA.fcsc ^-g^ 

Q* qb + ai/osc 

The Qesc/Q* can be called "effective" escape fraction and 
depends on / csc and T e through as. Note that the es- 
cape fraction / csc is the fraction of stellar and nebular 
LyC photons which escape from a galaxy relative to all 
(stellar+nebular) LyC photons produced in the galaxy (see 
eq. [7]). On the other hand, the "effective" escape fraction 
is the fraction of escaping stellar+nebular LyC photons rel- 
ative to the photons produced by only stars. 

Figure 8 shows the number fraction of LyC photons 
escaped relative to those produced by stars, or "effective" 
escape fraction, as a function of / eS c given by equation (8). 
The escaping number fraction is always larger than / esc but 



T e (10 4 K) 0.5 1 



a B /a A 0.669 0.626 0.573 0.502 
ai/a A 0.331 0.374 0.427 0.498 

(eneb) (eV) 16.41 16.61 16.98 19.12 



Table 3. Mean energies of stellar and escaping Lyman continua. 



Model 



BCD 



(e*> (eV) 22.66 23.00 26.85 32.18 

(ecsc) (eV) a 21.53 21.81 24.94 29.27 
(ecsc) (eV) 6 20.62 20.85 23.40 26.93 

a T = 10 4 K, / osc = 0.5 
b T c = 10 4 K, / esc = 0.1 



never exceeds unity. That is, the total number of escaping 
LyC photons is always smaller than the number of LyC pho- 
tons produced by stars. Nebular LyC compensates only a 
fraction of stellar LyC photons which absorbed by nebulae. 
Since the nebular LyC has a strong peak at the Lyman limit, 
however, the luminosity density near the limit is significantly 
enhanced as discussed in 53.2. 



3.4 Mean energy of escaping Lyman continuum 

Because of modification in spectra by nebular emissions, 
the mean energy of LyC photons escaping from galaxies is 
changed from that in the stellar spectrum. Let us define the 
mean energy of escaping LyC photons as 



(Cesc) 



-LyC 



L h J c + L 



LyC 
neb 



(9) 



Qcsc Q* ~\~ Qneb 

where L LyC s are frequency integrated LyC luminosities of 
escaping, stellar, or nebular radiations. The equal is true 
for frequency independent / osc . The Qs are already intro- 
duced in equations (1) and (7). If we define mean ener- 
gies of stellar and nebular LyC as (e») = L^ yC /Q* and 
(e ncb ) = L^/Qaeb, equation (9) is reduced to 

(e CS c) =/9(e.) + (l-/3)(e»eb>, (10) 
where 

p = («*) + (2L) fme . (ii) 

V qa / V a A ) 

Note that (e ncb ) = TncbV^i' where 7^ c y b c is the frequency 
integrated emission coefficient of nebular LyC. Therefore, 
the mean energy of escaping LyC is given by the internal 
division of mean energies of stellar and nebular LyC with 
weights P and 1 — /3, and thus, the escaping mean energy is 
always smaller than the stellar one if / eS c < 1. 

In Table 2 we summarize the coefficients for j3 in equa- 
tion (11) and the mean photon energy of nebular LyC as 
a function of T c . Table 3 gives a summary of mean photon 
energies of stellar LyC for models A-D in the first row. The 
second and third rows in Table 3 show mean photon energies 
of escaping stellar+nebular LyC for two cases of T c = 10 4 
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Figure 9. Spectra relative to the stellar one and comparison 
with results from Cloudy. The electron temperature To = 10 4 K 
and the age of 1 Myr after an instantaneous burst for the stellar 
population model D are assumed. The dotted line is the escaping 
stellar+nebular spectrum with a constant escape fraction / csc = 
0.5 which is the same as in Figures 3 and 4. The long-dashed 
line is the escaping stellar+nebular spectra for which we assume 
the neutral hydrogen optical depth at the Lyman limit tl =0.68 
and at shorter wavelengths to be oc A 3 . The solid and short- 
dashed lines are the cases calculated with Cloudy 08.00 (Ferland 
et al. 1998) for hydrogen only case and hydrogen and helium case, 
respectively. The Lyman limit optical depth for the hydrogen only 
case is matched to 0.68. 



K and / eS c = 0.5 or 0.1. We find that the energy reduction 
is about 5-15% for these cases and its maximum is 10-20% 
at the limit of / csc — » 0. 



4 DISCUSSION 

4.1 Wavelength dependent opacity and 
comparison with the code Cloudy 

As described in §2.2.1, LyC opacity through a clumpy ISM 
is independent of wavelength. As an opposite limit, we here 
consider a homogeneous ISM where the opacity keeps the 
wavelength dependence of the photo-ionization cross section 
of hydrogen: cx A 3 . 

Figure 9 shows escaping stellar+nebular spectra rela- 
tive to the stellar one in a homogeneous ISM (LyC opacity 
depends on the cube of wavelength) for the model D. The 
Lyman limit optical depth is set to be 0.68. This opacity 
corresponds to / csc = 0.5 at the Lyman limit and is applied 
to both stellar and nebular LyC (long-dashed line). If we 
compare the case in a clumpy ISM presented in the pre- 
vious sections (dotted line), we find that the Lyman limit 
bump predicted in a homogeneous ISM is weaker than that 
in a clumpy ISM. This is because higher energy LyC escapes 
easily and does not convert to nebular emission in the wave- 
length dependent case (i.e. homogeneous ISM). To make a 
strong bump, we need to convert higher energy LyC to the 
Lyman limit by clumpy nebulae. 

We can use the public photo-ionization code, Cloudy 
jFerland et al.| [l998) for a homogeneous ISM. To check our 
nebular calculations, we compare our results with the code 



Cloudy in Figure 9. We input the following condition to the 
code: a point source with the production rate of hydrogen 
LyC photons of lfj 52 01 s~\ which corresponds to 10 4 M 
in the stellar mass for the model D, resides at the center of 
spherical symmetric gas with the hydrogen atomic number 
density of 1 cm -3 and with the temperature of 10 4 K. We set 
the inner and outer boundaries to be 10 20 cm and 10 21 cm, 
respectively. The outer boundary is determined to achieve 
a Lyman limit optical depth of 0.68. This case is shown 
by the solid line in Figure 9. We find a good agreement 
between our calculation (long-dashed) and Cloudy (solid). 
However, Cloudy predicts somewhat stronger bump at the 
Lyman limit. This is because the optical depth for nebular 
LyC is smaller than that for stellar one since the nebula 
surrounds the source, whereas we set the same optical depth 
for both continua. 

The short-dashed line in Figure 9 is the case where the 
spherical gas consists of hydrogen and helium. Other set- 
tings are the same as the hydrogen only case. The resultant 
continuum is not very different from the hydrogen only case. 
Thus, the effect of helium is small for the shape of the SED. 
This point is also discussed in the next subsection. 



4.2 Effect of helium 

We can derive equations for helium nebular emissions like 
equation (4) for hydrogen. The strength of the helium emis- 
sions relative to that of hydrogen can be approximated to 



7? 



(12) 



where X is neutral helium (He I) or singly ionized he- 
lium (He II) and other quantities have the same mean- 
ings as in §2 but with superscript describing which atom 
or ion. In the four models of A-D considered here, we found 



0.4-0.7 and Q^ cll /Q™ 



4 x 10-0.1 with a 



larger value for a smaller metallicity l|Schaererll2002l . [20031 1 . 
The recombination rates of helium are Qr 6 ' w a™ and 

Hell 



Aot™ (|Osterbrock fc Ferlandl 120061). According to 



the em ission coefficients shown in lOsterbrock fe Ferland I 
(2006), we found that the helium contribution is negligible in 
the wavelength range of 800-2000 A, except for the model D 
where the He II bound-free Paschen continuum contributes 
about 20% around 2000 A. This smallness of the effect of he- 
lium on continuum can be confirmed in Figure 9 if we look 
at the result by Cloudy with hydrogen and helium (short- 
dashed line). However, the effect on recombination lines may 
not be negligible. 



4.3 Lyman "bump" as a tool to find exotic stellar 
populations 

As shown in Figure 4, we predict that galaxies dominated by 
very young (~ 1 Myr) and massive stellar populations can 
have a Lyman limit "bump" not "break" expected in the 
stellar spectra. The strength of the bump depends on the 
hardness of the LyC (see also Figure 6). If galaxies contain 
very massive (~ 100 M@) EMP stars or metal- free (Pop-Ill) 
stars, the Lyman limit bump becomes very strong. This fact 
provides a new method to find such exotic stellar popula- 
tions in galaxies at z < 4 although the IGM attenuation 
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hamp ers us to us e this method for galaxies at z > 4. In- 
deed, llwata et al.l (|2009Tj have found very interesting LAEs 
at z — 3.1 which emit extremely strong LyC. We will ap- 
ply the model developed in this paper to the LAEs in the 
forthcoming paper and discuss their nature in detail. 

4.4 Effect on reionization 

Let us examine if escaping nebular LyC affects the reioniza- 
tion process. The reionization history is regulated by the 
product of escape fraction and star formation efficiency, 
/esc/*, if we assume a cosmo logical structure formation sce- 
nario (|Wvithe fc Loebll2003l ). As the first effect of escap- 
ing nebular LyC, this / esc should be replaced by "effective" 
escape fraction, Q csc /Q*, which is introduced in equation 
(8) and Figure 8. From equations (8) and (11), we obtain 
Qesc/Q* ~ fcsc/P and find that the "effective" escape frac- 
tion is at most a factor of 1.6 larger than / osc for T c = 10 4 
K. As the second effect, the nebular LyC modifies the LyC 
spectral shape and reduces the mean photon energy of es- 
caping LyC. This leads to a decrease of the mean free path 
of LyC photons in the IGM. However, the reduction of the 
mean photon energy is at most 10-20% relative to that of 
the stellar LyC as discussed in §3.4 and in Tables 2 and 
3. Therefore, the nebular LyC escape itself may not give a 
significant impact on the reionization process. 



5 CONCLUSION 

The discovery of z — 3.1 LAEs emitting extreme ly strong 
LyC by the Subaru telescope l|lwata et al.l l2009h requires 
us to examine LyC emissivity and spectrum of galax- 
ies more carefully. As an attempt, this paper has exam- 
ined the effect of an escape of nebular LyC on SEDs 
of galaxies. The nebular LyC mainly emitted by bound- 
free transitions of hydrogen is usually assumed to be ab- 
sorb ed within nebulae, so-called "on-the-spot" approxima- 
tion (|Osterbrock fc Ferland 2006). However, we should con- 
sider its escape if st ellar LyC escapes from galaxies as 
observed at z ~ 3 JSteidel. Pettini. fc Adelbergerl l200ll ; 
IShaplev et"alll2006t llwata et al.ll2009h . 

Since the bound-free nebular emission has strong peaks 
at the limits of Lyman, Balmer, Paschen, etc. (see Figure 2), 
we expect a bump at the Lyman limit due to escaping neb- 
ular bound- free LyC (see Figure 3). This bump boosts lumi- 
nosity density at 900 A. As a result, 900 A-to-1500 A lumi- 
nosity density ratio can become larger than that expected 
in stellar SEDs. Indeed, the ratio has a maximum not when 
/esc = 1 but when / osc = 0.6 (see Figure 6). The strength of 
the bump depends on T c in nebulae (see Figure 5) , hardness 
of stellar LyC (i.e. metallicity, IMF and mass range, and age; 
see Figures 3 and 7), and IGM attenuation (see Figure 4). 
On the other hand, the total number of LyC photons which 
escape from galaxies is always less than that produced by 
stars (see Figure 8). In fact, the Lyman limit bump is a re- 
sult of the radiation energy re-distribution by nebulae. The 
energy re-distribution reduces the mean energy of escaping 
LyC photons, but it is at most 10-20% reduction relative to 
the mean energy of stellar LyC (see Table 3). 

We have assumed wavelength independent / osc (or LyC 
opacity through the ISM) to obtain our main results. This 



should be suitable for a clumpy ISM (see schematic Figure 1 
and discussion in §2.2.1). We have also calculated a wave- 
length dependent case which corresponds to a homogeneous 
ISM, and have found that the wavelength dependent case 
predicts a weaker Lyman limit bump (see Figure 9). In ad- 
dition, we hav e compared the case with the photo-ionization 
code Cloudy (jFerland et al.lll998l ), and have found a good 
agreement (see Figure 9). Moreover, we have omitted he- 
lium, but it is justified by a discussion in §4.2 and by a 
comparison with Cloudy (see Figure 9). 

The most interesting implication from our results would 
be the prediction of Lyman "bump" galaxies (see Figure 4 
and §4.3); galaxies containing very young (~ 1 Myr), very 
massive (~ 100 Mq), and extremely metal-poor (or metal- 
free) stellar populations show a prominent bump just below 
the Lyman limit. This would be a new indicator to find 
such exotic stellar populations. Although IGM attenuation 
hampers us to use this method for z > 5 galaxies, we can 
use it for z < 4 galaxies. 

Finally, we have discussed effects on the reionization 
briefly (§4.4). Escaping nebular LyC has two effects: acting 
as an additional LyC source and reducing the mean energy 
of escaping LyC photons. The first effect increases the escape 
fraction effectively, but the enhancement is at most a factor 
of 1.6 (for T c = 10 4 K) and the LyC photon escape rate never 
exceeds the production rate by stars (see Figure 8). The 
second effect decreases the mean free path of LyC photons 
in the IGM. However, the reduction of the mean photon 
energy is at most 10-20% relative to the stellar LyC (see 
Table 3). Therefore, the nebular effect on the reionization 
may be small. 
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